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PREFACE 


This  study  involved  tasks  related  to  three  WES  programs. 
The  initial  program,  DASA  Subtask  14.020  Deep  Water 
Wave  Study  is  presented  in  Volume  I of  this  report.  The 
latter  programs  DASA  Subtasks  14.080  Mono  Lake  Field 
Experiment,  and  14.096  Shallow  Water  Wave  Study  are 
covered  in  Volumes  II  and  III  respectively. 
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2.  THEORETICAL  DEVELOPMENTS  RELATED  TO  WATER  WAVES 
GENERATED  BY  EXPLOSIONS  IN  DEEP  WATER 

This  section  has  been  written  with  the  objective  of  tracing  the 
historical  development  of  the  multitude  of  papers  which  have  contributed 
to  the  theory  presently  being  applied  for  predicting  water  waves  generated 
by  deep  water  explosions  in  a constant  depth  environment.  A general 
theoretical  presentation  is  included  and  various  special  cases  are.  presented 
with  appropriate  reference  to  original  contributions.  The  mathematical 
development  is  carried  out  in  detail  to  minimize  the  difficulties  for  readers 
not  thoroughly  familiar  with  the  theory.  A discussion  of  the  initial  conditions 
required  for  performing  predictions  is  given  along  with  an  evaluation  of  the 
parameters  using  the  experimental  data  resulting  from  this  program.  Finally, 
possible  modifications,  extensions,  and  applications  of  the  theory  are 
summarized  in  a detailed  format  which  can  be  used  as  a guideline  for 
further  theoretical  development. 


The  linear  theory  of  water  waves  generated  by  explosions  in  deep 
water  of  constant  depth  consists  of  solving  the  Laplace  equation  in  the 
appropriate  coordinate  system  subject  to  the  linearized  free  surface 
condition,  the  bottom  condition  and  initial  conditions,  equivalent  to  the 
position  and  normal  velocity  or  velocity  potential  of  a fluid  surface  at 
one  boundary.  Several  sets  of  initial  conditions  can  be  considered: 

1.  An  initial  deformation  and  velocity  field  on  the  free  surface 

2.  An  initial  shape  of  the  free  surface  and  an  applied  pressure 

3.  An  initial  deformation  and  velocity  field  on  the  bottom. 
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The  first  theoretical  contributions  on  the  theory  of  waves  pro- 
duced by  a local  disturbance  in  deep  water  can  be  found  in  the  classical 
memoirs  of  Cauchy  (Ref.  1)  and  Poisson  (Ref.  2).  These  works  were 
performed  independently  and  yield  essentially  equivalent  results  for  two- 
dimensional  problems.  A mathematical  analysis  of  both  papers  was 
performed  by  Burkhardt  (Ref.  3),  and  Lamb  (Ref.  4)  has  reviewed  the 
theory  and  presented  it  in  a consistent  manner  including  contributions 
of  his  own.  The  two  principle  shortcomings  of  these  early  papers  were 
1)  the  approximations  necessary  for  evaluating  the  integrals,  and  2) 
confining  the  initial  elevation  or  impulse  to  the  vicinity  of  the  origin. 

Lamb  (Ref.  4)  improved  the  mathematical  approximations  by  observing 
that  he  could  make  use  of  a tabulation  of  the  power  series  used  by  Rayleigh 
(Ref.  5)  which  appeared  in  the  theory  of  Frenels  diffraction  integrals. 

The  mathematical  difficulties  which  burdened  the  early  works 
were  removed  by  Kelvin  (Ref.  6)  when  he  developed  the  method  of  stationary 
phase  for  approximating  the  integral  at  large  distances  from  the  source  of 
the  disturbance.  Many  authors  made  contributions  to  the  understanding  of 
these  early  classical  papers  with  the  development  of  the  concept  of  group 
and  phase  velocities,  a description  of  the  phenomena  of  waves  in  a dis- 
persive medium,  and  various  mathematical  methods  of  approximation. 

A few  of  the  more  notable  of  these  developments  were  due  to  Watson  (Ref.  7), 
Burnside  (Ref.  8),  Dawson  (Ref.  9).  Terazawa  (Ref.  10),  Pidduck  (Ref.  11) 
Thompson  (Ref.  12),  PidducKRef.  13),  Stokes  (Ref.  14),  and  Havelock 
(Ref.  15). 

The  early  classifical  results  above  were  supplemented  by  several 
papers  in  Japan  on  the  waves  produced  by  a bottom  deformation,  the  impetus 
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for  the  consideration  of  this  problem  being  the  phenomena  of  tsunami 
waves  produced  by  seismic  disturbances  of  the  sea  bed.  The  earlier 
papers  were  by  Hasegawa  and  Sano  (Ref.  1 6) , Sano  (Ref.  17),  and 
Terazawa  (Ref.  18). 

Further  theoretical  developments  relating  to  explosion  waves 
were  practically  non-existent,  in  England,  Japan,  and  the  U.S.  , until 
World  War  II  when  a need  for  more  detailed  knowledge  about  underwater 
explosion  characteristics  precipitated  a rapid  development  of  the  theory. 

At  this  time,  the  work  in  Japan  was  still  related  to  the  tsunami  problem, 
some  of  the  more  notable  contributions  being  made  by  Nomura  (Ref.  19), 
Segawa  and  Kanai  (Ref.  20),  Takahashi  (Ref.  21,  22,  23),  Ichiye  (Ref.  24), 
Honda  and  Nakamura  (Ref.  25),  Nakamura  (Ref.  26),  Unoki  and  Nakano 
(Ref.  27,  28,  29),  and  Kajiura  (Ref.  30).  This  latter  paper  by  Kajiura  is 
the  most  notable  in  that  it  presents  the  entire  linear  theory  in  a self- 
contained  consistent  manner  which  can  be  applied  for  sea  bed  disturbances, 
surface  disturbances,  pressure  distributions,  and  an  initial  velocity 
distribution.  A general  approach  is  taken  to  the  problem  and  it  is  solved 
by  the  use  of  a time  dependent  Green's  function  with  several  special  cases 
being  worked  out  for  the  asymptotic  solution.  This  paper  and  related 
references  will  be  referred  to  later  during  the  development  of  the  mathe- 
matical formulation  and  solution. 

Extensive  work  was  performed  both  in  England  and  the  United  States 
during  and  immediately  after  World  War  II  on  what  is  commonly  referred 
to  as  the  wave  "generation  mechanism".  This  phrase  means  the  early 
explosion  behavior,  i.e.  , the  propagation  of  the  detonation  wave,  production 
of  explosive  gases,  formation  of  an  underwater  bubble,  and  bubble  oscillations 
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and  migration.  Many  of  these  papers  are  collected  in  Refs.  31  and  32 
and  much  of  the  work  also  appears  in  Cole's  well  known  book  (Ref.  33). 

In  addition,  Snay  (Refs.  34-36)  has  made  notable  contributions  over  the 
past  15  years  to  the  understanding  of  explosion  bubble  phenomena.  Bubble 
phenomena  will  not  be  discussed  in  this  report  due  to  the  vast  amount  of 
literature  on  the  subject  and  the  fact  that  state  of  the  art  prediction  tech- 
niques do  not  require  any  knowledge  of  bubble  phenomena. 

With  the  advent  of  nuclear  explosions,  a new  and  perhaps  catas- 
trophic problem  was  foreseen:  the  possible  generation  of  a huge  tsunami- 

like surface  wave  which  might  produce  extensive  coastal  damage  along  the 
entire  East,  West  or  Gulf  Coasts  of  the  United  States  from  a single  detonation 
in  deep  water.  This  problem  was  attacked  through  a comprehensive  pro- 
gram formulated  and  monitored  through  the  Department  of  Defense.  Although 
mathematical  models  had  been  developed,  the  applicability  of  these  models 
was  unknown  as  well  as  a realistic  appraisal  of  the  initial  conditions. 

Stoker  (Ref.  37)  gave  the  Green's  function  for  a three-dimensional  dis- 
turbance in  water  of  constant  depth.  The  first  successful  attempt  to 
formulate  a realistic  theory  applicable  to  the  prediction  of  explosion  waves 
was  published  by  Kranzer  and  Keller  (Ref.  38).  Prins  (Ref.  39)  performed 
a laboratory  experiment  in  a two-dimensional  tank  (i.e,  no  radial  spreading) 
where  he  elevated  or  depressed  the  water  surface  at  one  end  of  the  tank 
and  measured  the  subsequent  wave  profile  at  several  stations  for  various 
ratios  of  water  elevation  to  water  depth.  The  initial  major  efforts  of  the 
DOD  wave  program  consisted  of  a review  of  available  theory  and  data, 
both  HE  and  nuclear,  and  two  experimental  test  series.  The  review  was 
performed  by  Kaplan  (Ref.  40)  of  Broadview  Research  Corporation  (later 
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to  become  United  Research  Services).  Many  shortcomings  were  evident 
from  this  study,  especially  with  respect  to  well-controlled  experiments 
for  obtaining  data  on  both  the  wave  generation  mechanism  (with  the  objective 
of  relating  initial  conditions  of  the  theoretical  solution  to  the  actual  mode 
of  generation)  and  the  wave  train  produced  by  deep  water  explosions  in 
water  of  constant  depth.  The  experimental  test  series  on  the  wave 
generation  mechanism  was  performed  at  a small  scale  under  reduced 
atmospheric  pressure  in  a tank  filled  with  oil.  Hemispherical  charges 
were  detonated  at  the  glass  wall  of  the  tank  and  high  speed  motion  pictures 
were  taken  of  the  subsequent  bubble  and  plume  formation.  The  tank  was 
not  extensive  enough  for  measurements  of  the  wave  profile  to  be  taken 
at  any  significant  distance  from  the  charge.  However,  a very  good  quali- 
tative understanding  of  the  method  of  wave  generation  can  be  obtained 
from  a careful  observation  of  these  films.  This  test  series  is  reported 
on  in  a series  of  eight  reports  by  Kaplan,  et  al.  (Refs.  41-48).  The  other 
test  series  was  performed  by  the  Waterways  Experiment  Station  and  this 
volume  constitutes  the  final  report  o f that  test  series.  Earlier  reports 
of  the  data  from  the  various  field  tests  and  an  analysis  of  propagation 
characteristics  may  be  found  in  Refs.  49-52.  These  field  tests  consisted 
of  detonating  1/2  lb  to  1000  lb  charges  of  TNT  and  measuring  the  water 
waves  at  many  stations  in  water  of  constant  depth.  Due  to  the  compre- 
hensive planning  of  these  tests,  the  data  has  resulted  in  an  evaluation  of 
the  initial  conditions  for  the  mathematical  model  as  a function  of  yield 
and  submergence  depth  resulting  in  predictions  of  a known  reliability 
for  HE  charges  over  many  orders  of  magnitude  in  charge  weight.  Extra- 
polations to  nuclear  explosions  are  still  more  uncertain;  however,  it  is 
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generally  accepted  that  periods  can  be  predicted  within  20%  and  ampli- 
tudes at  least  within  a factor  of  2 and  probably  considerably  better.  The 
uncertainty  stems  from  a serious  lack  of  deep  water  nuclear  tests  and 
the  absence  of  reliable  wave  measurements  from  the  few  that  were  per- 
formed. The  present  state  of  the  art  prediction  technique  for  deep  water 
wave  trains  and  the  answer  to  associated  questions  relating  to  analytical 
approximations  can  be  found  in  a series  of  reports  by  Whalin.  et  al,  from 
Interstate  Electronics  Corporation  (Ref.  53-59).  A discussion  of  the 
applicability  of  the  theoretical  methods  to  areas  near  the  source  of  dis- 
turbance is  discussed  by  Whalin  (Ref.  60).  A relatively  early  exposition 
on  impulsively  generated  water  waves  was  authored  by  Van  Dorn  (Ref.  61) 
and  some  measurements  from  nuclear  detonations  are  discussed  in  Ref. 

62  and  63.  Van  Dorn  also  reported  on  a field  program  involving  14,  400 
lb  HBX  charges  and  the  relationship  between  the  measured  profiles  and 
those  that  could  he  reproduced  through  theoretical  considerations  (Ref.  64). 

It  should  be  noted  that  with  the  analysis  of  data  from  this  test  series 
and  the  availability  of  data  from  the  1965  and  1966  Mono  Lake  tests,  the 
state  of  the  art  in  predicting  deep  water  explosion-generated  waves  is 
undergoing  very  rapid  change  and  modification.  A very  interesting  project 
which  provided  for  an  assessment  of  the  state  of  the  art  Decame  available 
through  the  prediction  by  Le  Mehaute  and  Whalin  (Ref.  65)  of  the  waves  to 
be  produced  from  several  explosions  in  Mono  Lake.  Preliminary  com- 
parisons with  the  subsequent  experimental  data  indicate  an  approximate 
accuracy  of  20%  in  amplitude  and  10%  in  period.  Another  report  by  Whalin 
(Ref.  66)  on  this  same  project  summarized  the  theoretical  background  for 
the  predictions  and  indicated  several  modifications  and  extensions  to  the 
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theory.  The  following  section  of  this  report  fully  contains  those  findings 
of  Ref.  66  and  considerable  supplementary  information.  An  interesting 
application  of  the  state  of  the  art  as  of  one  year  ago  can  be  found  in 
Refs.  67-69  where  LeMehaute,  Whalin,  et  al.  predicted  coastal  damage 
from  waves  near  Los  Angeles,  Providence  and  New  Orleans  from  offshore 
underwater  and  near  shore  air  detonations. 

The  majority  of  research  being  conducted  at  this  time  is  devoted 
to  a better  understanding  of. the  modification  of  a wave  train  as  it  propa- 
gates over  the  continental  slope  and  shelf,  and  significance  of  three- 
dimensional  effects  on  run-up,  the  run-up  from  nonbreaking  waves,  shelf 
and  harbor  resonance,  structural  damage,  and  various  other  applications 
of  the  theory  to  specialized  problems.  The  primary  need  for  improving 
deep  water  predictions  is  a well  documented  nuclear  test.  Various 
modifications  and  extensions  also  should  be  carried  out  over  the  next 
several  years  as  indicated  in  the  following  sections  of  this  chapter. 
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The  theoretical  formulation  of  the  problem  will  be  presented  in 
a general  form  as  given  by  Kajiura  (Ref.  30)  and  Whalin  (Ref.  66);  however, 
most  details  of  the  mathematical  development  are  presented  as  well  as  a 
discussion  of  the  contribution  of  various  terms  and  the  physical  meaning 
of  assumptions  in  both  the  original  formulation  and  the  method  of  solution. 
The  energy  of  the  source  of  disturbance  also  is  given  in  detail. 

A multitude  of  classical  and  recent  theoretical  papers  was  referenced 
in  the  preceding  section.  Regardless  of  their  particular  objectives,  one 
facet  is  common  to  all  of  these  papers,  namely  the  use  of  a linear  theory 
for  surface  waves  in  an  inviscid  and  incompressible  fluid.  The  assumption 
of  constant  depth  is  almost  a necessity  for  obtaining  closed  form  analytical 
solutions  and  also  is  common  to  the  majority  of  publications  related  to 
this  problem.  The  principle  difference  in  these  references  consists  of 
various  assumptions  for  initial  conditions,  the  use  of  two-  or  three- 
dimensional  coordinate  systems,  and  different  methods  of  approximation 
for  obtaining  an  analytical  solution.  The  initial  formulation  and  solution 
will  include  all  previous  publications  as  special  cases  and  the  solutions 
for  special  cases  will  be  discussed  in  detail. 

Our  assumptions  will  consist  of  the  following: 

1.  The  water  is  incompressible 

2.  The  fluid  is  non-viscous 

3.  The  flow  is  ir rotational 

4.  The  fluid  domain  is  of  constant  depth,  d,  and  infinite  in  extent 

5.  The  first  order  of  approximation  is  assumed  for  the  kinematic 
and  dynamic  conditions  at  the  free  surface. 
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The  origin  of  the  Cartesian  coordinate  system  (x',  y',  z')  is  at 
the  undisturbed  free  surface  with  the  vertical  axis  z’  being  taken  positive 
upwards.  All  variables  are  transformed  to  dimensionless  form  (x=x'/d, 
y = y'/d,  z = z'/d,  t = t'  Jg/d,  T)  = t|'/d,  cp  = cp'/d  J gd  , p = (p'/p)/gd,  etc). 

The  kinematic  and  dynamic  conditions  at  the  free  surface  are 


cp  -n  -cp  +cp  -n  +TT.  = 0 
x 'x  z y y x 

A .1,2.  2 ^ 2,  x n 
T1  + CPt  + 2 + 'Py  + 'Pz  ) + P = 0 

With  assumption  (5)  on  the  previous  page,  these  conditions  are  linearized 
and  become 


<pz  = V 2 = 0 

CPt  = -Tl-p,  z = 0 


The  bottom  condition  is 


^n  = WB’  z = - h (x,  y) 


which,  from  assumption  (4),  becomes 


CD  = W 


z B 


z'=-d 


or  in  dimensionless  form 


rPz=vB*  I = '1 

where  w'  is  the  assumed  bottom  velocity  corresponding  to  a bottom 
B 

deformation  and  cp’  is  the  velocity  potential,  which,  from  assumptions 
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(1)«  (2),  (3)  and  (4)  assures  the  existence  of  a single  valued  function 
cp  (x,  y,  z,  t)  from  which  the  velocity  field  can  be  derived  by  taking  the 
gradient 

v = grad  cp 

Since  the  velocity  components  are  given  by  the  above  equation,  it  follows 
directly  from  the  continuity  equation 

div  v = 0 

that  the  velocity  potential,  cp,  is  a solution  of  the  Laplace  equation 

H ^ =0 
d x d y 3 z 

and  therefore  cp  is  a harmonic  function. 

It  should  be  noted  that  linearization  places  certain  restrictions 
on  the  surface  and  bottom  deformations  that  may  be  assumed.  Specifically, 
both  should  be  small  compared  with  the  wavelength  X.'  and  water  depth  d. 
The  condition  z'\,2/d3  ~ 1 should  also  be  satisfied.  However,  we  will 
see  in  the  next  sections  that  for  practical  purposes  of  prediction,  the  defor- 
mation may  be  much  larger  than  one  would  normally  expect  for  the  theory 
to  be  valid. 

This  problem  will  be  approached  through  the  use  of  a time 
dependent  Green's  function,  G,  which  is  a solution  of  the  Laplace  equation 

V2  G = 0,  0 > z > - 1 , tST 
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satisfying  the  free  surface  condition 


Gtt  + Gz  = °’  z = 0 


and  the  bottom  condition 


G =0,  z = - 1 
z 


The  mathematical  conditions  on  the  Green's  function  which  must  be 
satisfied  to  assure  the  existence  of  a unique  solution  are  that  G,  G^, 
G^,  Gz>  G^,  G^,  G^y,  and  Gtz  must  be  uniformly  bounded  for  every 
t at  x,  y — > oo  and  (G  - 1/R)  must  be  bounded  at  any  source  point 
(xQ.  yo.  zo)  where 


R2  = (x-xo)2  +(y-yQ)2  +(z-zQ)2 


At  time  t = T the  Green's  function  is  assumed  to  satisfy  the  following 
initial  conditions 


G = Gt  = 0,  z = 0 


The  above  equation  and  conditions  are  sufficient  to  determine  the  Green's 
function  uniquely.  For  a derivation  of  the  Green's  function,  in  the  case 
of  three-dimensional  motion  in  water  of  finite  depth,  the  reader  is  referred 


to  Stoker  (Ref.  37).  From  Stoker 


G<v  WTh y<«4>  * | jsarr  [sinha  0 -Is -si}  - sinh  0 {1+(z+s>} 

+ |l  - cos  y (t-T)}  co8g--  cosh  a (1  +z)  cosh  a (1  +zQ)]  da,  0 > z,  zq  > - 1 
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Upon  making  use  of  Green's  formula,  the  above  Green's  function,, 
and  symmetry  properties  of  G,  one  obtains  an  expression  for  cp^.  (x,  y,  z,  t). 
The  integration  of  cp^  from  0 to  t and  subsequent  substitution  of  cp 
into  the  boundary  conditions  gives  the  following  expression  for  the  free 
surface  elevation,  T] 


^+P  =T?ll[f¥z  - ^t^T-O  + [JpGTTtdT  +(pGTT)T=tj}z=Zn=0 


b 

+ [J  S",  rdT]z  =-l,  z=o]dS 


The  details  of  the  above  derivation  may  be  found  in  either  Stoker  (Ref.  37), 

Kajiura  (Ref.  30)  or  Whalin  (Ref.  66).  The  above  expression  for  r\  represents 
the  classical  solution  derived  in  a considerably  more  general  form. 

Lets  now  analyze  the  contribution  of  each  term  in  the  above  equation 

1.  The  first  term  represents  the  contribution  of  the  initial  velocity  j 

distribution,  cp  , assumed  on  the  free  surface. 
zo 

2.  The  second  term  represents  the  contribution  of  the  initial  defor- 
mation, *n JT _0 , of  the  free  surface. 

3.  The  third  and  fourth  terms  represent  the  contribution  of  the 

applied  pressure,  p. 

4.  The  latter  term  represents  the  contribution  of  the  bottom  deformation, 
w being  the  vertical  velocity  of  the  bottom  deformation. 
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2.3  Asymptotic  Solution  for  Special  Cases 


The  difficulties  involved  with  the  solution  of  the  above  equation 
depend  upon  the  area  where  one  wishes  to  obtain  an  estimate  of  the  wave 
profile  and  the  functional  dependence  of  the  initial  conditions  on  the  space 
and  time  coordinates.  The  solution  can  be  integrated  numerically  regard- 
less of  the  complexity  of  the  initial  conditions.  However,  in  the  case  of 
waves  generated  by  a single  explosion,  the  problem  is  axially  symmetric. 
This  results  in  the  space  integration  becoming  the  zero  order  Hankel 
transform  of  the  initial  mode  of  disturbance  or  the  integral  of  the  dis- 
turbance from  0 to  t.  If  one  cannot  perform  the  Hankel  transform  ana- 
lytically, the  approximate  integration  method  developed  by  Whalin  (Ref.  60) 
may  be  used  to  accurately  compute  the  amplitude  time  history.  A closed 
form  solution  may  be  obtained  by  use  of  the  method  of  stationary  phase  to 
approximate  the  integral  over  the  wave  number.  The  following  examples 
illustrate  the  closed  form  solution  for  several  cases  of  the  various  types 
of  initial  conditions. 

Upon  computation  of  the  appropriate  partial  derivatives  and 
expressing  the  solution  in  cylindrical  coordinates,  one  obtains  the  following 
relationships . 


Gt|z=0,  zo=0,  T=0 
GtT | Z = 0 , ZQ  = 0> 


00 


-L2-  j ( o r)  sin  Y t do 
Y o 


oo 


j 2 0 Jq  ( o r)  cos  Y t d 0 


JL  J 


CD 


*TTt  | : 


5=0,  z =0  = ■ J 2 0 Y J0  ( 0.1)  sinv  (t-T) 


d o 


CO 


TT  Z 


n o , = 2o  J (Or)d(J 
=0,  z =0,  t=t  J o'  ' 


oo 


! 


cti, 


= -i , z=o  = J Jo(ar)  sinY  (t’T)da 


oo 


JQ  (a  r)  = Jo  ( a r)  Jq  ( a .rQ)  + 2 £ Jn  ( 0 r)  Jn  ( ° ro)  COS  n " 0o) 

n=l 

Note  that  upon  integration  with  respect  to  the  source  variable  0Q  (assuming 

axial  symmetry  for  all  initial  disturbances) 

2tt 

J V o7)dS,o  = 2"Jo,0r)  Jo(oro> 
o 

oo  oo 

'l  + P - J [Jf  J0Cor)J0(or0)»mYtdo  cp  _q 
O o 'o’ 

t 

+ j J0(ar)  Jo  (a  ro)  cos  yt  d a r\  (ro)jT=0 
o 

t oo 

+ f p (r  , t)  f - a Y j ( a r)  J ( a r ) sin  Y (t  - T)  d o d t 
j O O Q 

o ° 

oo 

+ (P  (ro’  T)  I a‘  JQ  (ct  r)  Jo  ( ct  ro)  d a)r=t 
o 

t oo  ' 

- j J7^h-si,,Y(t-T)d0  =.1dTJr»dr° 

o O O z A 


Note:  If  the  initial  surface  deformation  is  defined  in  general  by  T)o(ro>  t) 
then 


% (ro’  T)1  T =0  = % (ro> 
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Vk-o-r-o  = aT  'To  <V  T)ko  ’ % <V°> 

O O T 


Now  consider  the  following  special  cases: 


Case  1:  The  applied  pressure  at  the  surface  is  impulsive  at 


T = 0 


P (r0>  T)  = p (ro)  6 (0  ) 


therefore 


j P (ro,  x)  GTn  dT  + p (rQ,  r)T=t  = ID  (rQ) 


All  other  terms  are  identically  zero  of  the  equation  since  = 0,  T|  (rQ)  = 0, 

o 

Wg  = 0.  Therefore 

oo  oo 

H(r.t)  = j a Y JQ  (a  r)  sin  V t j IQ  (rQ)  JQ  ( o rQ)  rodrQdo 


Upon  inspection  one  notices  that 


f I (r  \ j ( CT  r ) r dr  =1  (a)  = zero  order  Hankel  transform 

J o ' o o ' o o o o 


of  the  initial  impulse. 
Therefore 


n (r,  t)  = j o Yl0(o)  JD  (a  r)  sin  Y t do 


The  above  integral  can  be  evaluated  by  the  method  of  stationary 
phase  subject  to  the  following  assumptions 
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1)  The  frequency  of  oscillation  of  IQ(a)  is  low  compared  with 

that  of  J (or)  and  sin  vt. 

o ' 

2)  JQ(  a r)  may  be  approximated  by  the  first  term  of  its  asymptotic 
expansion  for  large  argument.  Note:  Additional  terms  in  the 
asymptotic  expansion  will  lead  to  a slightly  better  approximation; 
however,  the  increase  in  accuracy  is  small  compared  with  the 
uncertainty  in  choosing  the  initial  conditions. 

We  now  proceed  with  the  evaluation  of  the  integral,  preceded  by  the  following 
theorem. 


Theorem:  If  an  integral  is  of  the  form 


. X(*> 


T f (z) 


where  X (z)  is  of  bounded  variation,  T is  a large  positive  constant,  Tf(z) 
is  analytic  in  (A,  B),  and  f (z)  = cp  + if  , where  cp  = 0 and  f = 0 has  a 


unique  solution  xq  in  (A,  B),  then 


j X (z)  eTf<z)  d!  » x<*0> 


i(T*<*0)  +J»qn*"*0}+  L 


V T 'll"  (*„) 


Proceeding  with  the  evaluation,  we  make  the  following  assumptions 
1)  OY  r (a  ) has  a low  frequency  of  oscillation  relative  to 


JQ  ( o r)  and  sin  y t. 


2)  J (cr  r)  is  approximated  by  the  first  term  in  the  asymptotic 


V"1 


expansion  for  large  argument: 


hf-r  -s(ar-J)  = 


2 e 


i ( c r - S’)  + e-i(o  r - J) 


TT  CT  r 


Upon  substitution  one  obtains 


r'(r*  * I 7 aYIo(o } 


i(o  r +yt  - j)  i(o  r -Yt  - J) 


-i(o-r-yt-J)  -i(or+Yt-J) 


+ 


The  expression  for  r|  (r,  t)  is  row  in  a form  which  can  be  evaluated 
by  the  method  of  stationary  phase  where  cp  = 0,  ♦ = ^ (or  - Yt  - tt/4) 
and  the  stationary  phase  points  crQ  are  solutions  of  i) i' (a  ) = 0,  where 

2 

f(a)  = + r + ^nh  P + a sech  p t 

Vatanho 


2 ,2 


\ _ n Ttanh  cT  + a sech  al  , sech  a - a sech  a . tanh  cr.  \ 

t ( a ) - 17  ; 572 — *—  + T7? 1 t 

(otanh  a)  ' (a  tanh  cx.)i/2  J 


The  stationary  phase  points  are  defined  by 


cp(a  ) = 2 


_ 1 / / tanh  Q 


cosh  o 


2 v tanh  o J t 


Only  the  second  and  third  integrals  have  real 
stationary  phase  points.  In  the  second  integral 


1(a)  - a r-vt-^- 


V(a)  = J* 


1 4/ 

tanh  o + q sech  q 
2 ,/"a  tanh  0 


...  . f 1 r tanh  n + a sech2  - ]2  , sech4  0 - a sech4  o tanh  a \ 

[ o tanh  o J 


The  stationary  phase  points  are  defined  by 


1 / /tanh  g 1 J i 

®(a)-rw  2 cosh2  a Vtanh 


and  cp"(o  ) i8  positive. 
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In  the  third  integral 


Uo)  = -(ar  -Yt  -4) 
V (cr  ) ~ ~ { 


tanh  0 + a sech  a , 1 
- 1 r - ' t j 

2 V 0 tanh a 


H"(a>  = 


r 1 r tanh  a + a sech  CT-]“  , sech  0-0  sech  a tanh  0 1 
^ 4 [a  tanh  a]372  [atanha]1/2 


» 


The  stationary  phase  points  are  defined  by  the  same  functional  relationship 
as  in  the  previous  case,  however  iji"  ( 0 ) is  negative  in  this  case. 

Upon  application  of  the  previous  theorem 


i{{  ar  -yt-  J+J(+l)}  -if' a r - yt -J- J(-l)} 

+ e 


TT  TT 


^(r*  t)  = a Y IQ  (a  ) 7^7 


2 tt  e 


(o  )t 


4 i 


where  I'll"  (a  )|  has  been  replaced  by  (-cp'(o)  t).  Upon  substitution  of 
t = r/c o (a  ) and  simplification,  one  obtains 


2 I (a  )a  Y f- 

cp  ( 0 ) / 

ei(a  r -Yt)  +e-i(CT-r-Yt)^ 

V- 

2 IQ  (a  ) a Y 

0 cp1  (a  ) l 

4 i J 

[-  sin  (a  r - Yt)  j 

/ T(o)— 

r 'l 

- 0 co'  (0  ) 

Io(0)0. 


Jy(aTJanhZ  sin  ( 0 r . J Q tanh~  t) 


> 
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The  above  equation  reduces  to  Kranzer  and  Keller's  result  when  put  in  the 


form 


ri  (r,  t)  = 


Xo(0)  ° /co ( a ) tanh  a 


7 ) tanh  a . _ ,t  r, 

- ep-(— r 8in2n(f-\) 


where 


cp(  o)  = 


1 I I tanh  o 


cosh  0 


2 _ V tanh  o 


i0<°) 


= f I (r  ) J ( 0 r ) r d r 
J o ' o o o o o 


o tanh  a 


All  variables  above  are  in  dimensionless  form  as  previously  noted. 


Case  2:  The  applied  pressure  is  a function  of  time  from  0 to  t 


and  zero  thereafter 


p(ro.T),  T*T 


0 , T > T 


therefore 


j P <ro’  T)  GTTt 


dr  + [p(r0.  r)]T=t  = Jp<V  t)  GTTt  d 


oo 

= [ -0YJo(°r)J0(«.r0)  jp(ro.T)sinv(t-T)dTdO 


J " ~ T “o 
o 


i 


Let 


T* 

j P (rQ>  t)  sin  v (t  - t)  dT  = P (r  , t) 
o 

then 

oo  00 

Tl(r.t)  = - j aYJ0(a.r)  JP(r0,t)J0(ar0)r0dr0da 
o o 

Let 

00 

I P (r  , t)  J (a  r ) r dr  =P(o,t)=  zero  order  Hankel  transform  of  P (r  , t) 
J ' o o ' o o o ' ' ' o 

o 

CD 

T|  (r,  t)  = - f a Y P(a  ,t)  Jo  ( 0 r)  d a 
o 

The  above  integral  also  may  be  integrated  by  the  method  of  stationary 
phase  once  P(a,t)  is  known  and  provided  that  stationary  phase  points 
exist.  Otherwise,  if  one  cannot  integrate  directly  (which  is  unlikely)  then 
another  approximate  or  numerical  integration  technique  must  be  used. 


Case  3:  The  applied  pressure  is  non-zero  for  all  time  t 


P = P t) 


therefore 


oo  t 


1 r dr  - P(r,  t) 


^ = * li  P (ro’  GTTt  d T + LP  {l*0’  T)  GTT  JT=t  / rQ  d r0 


o o 


oo 


00 


00 

i r 


- J „J0l,r|  [-  I P (rD.  <V  *>  L < •%»  ro  d-oJ  O" 

W O L W W Q 

O O 

- P(r,  t) 


where 


P (ro»  = J P (rQ.  t)  sin  y (t  -T)dr 


(a.i)  = JP(r0,t)  Jo(aro)rodrc 
o 
OO 

( a,  t)  = p (r  , t)  J (ar  )r  dr 
J e ' o ' o ' o'  o o 


y-Kj 

(r,  t)  = j o Y Jq  (a  r)  j - P (a  , t)  + P d a - P (r,  t) 


Similar  comments  are  applicable  as  in  the  previous  case  with 
regard  to  the  integration. 

Case  4:  There  is  an  initial  surface  deformation  which  is  stationary 

at  time  zero.  Therefore 


r\  (r,  t)  = 


oo  oo 
p p 


j j a Jc  (a  r)  Jo  (oro)  cosyt  ^ (rQ)  d a rQ  d r£ 


oo  “ 

= I a jn  ( o.r)  COS  Yt  j r)Q  (rQ)  JQ  (a  rQ)  rQ  d 


r d a 
o 


; ^ ( a)  a JQ  (a  r)  cosyt  do 

J O 


Upon  evaluation  by  the  method  of  stationary  phase 
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Case  6:  An  initial  surface  deformation  and  an  initial  velocity  field. 

The  solution  in  this  case  is  merely  the  sum  of  the  solutions  for  Cases 


Case  7:  A bottom  deformation  of  constant  velocity  for  a 

finite  time  T*.  Therefore 


<Gt°Z  >T=0  = 


f 2 a sin  yt  d o 

WB  J Y cosh  cr 
o 


* 

T oo 


j CtT  dT 
o 
o 


ri  (or) 

, 1 -2 — — * 2 a[-  cosy(t-T)  do. 

B J J cosh  a L 


^ oJn(ar)Tr 

- _ w ! 2 i cos  Y(t  - t)  d r d a 

B i cosh  a J 

o o 


T q Jq  ( a r)  f - sin  Y (t  - t)1T 

= " 2 wb  J cosh~a  L Y J0 
o 


! da 


I*  ° [sin  y (t  - T*)  - sin  Yt]  d o 

I > > e h rV  *■ 


“ 2 wb  J v cosh  a 


Therefore 


“H  (r  > t)  - 2 Wg  j y cosh  a 


>)  - sinVt]  j ■•odrod':, 


“ n(o)  a . T»,  . si„Yt]  da 

= 2 wb  <1  Y cosh  a 
o 
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2. 4 Evaluation  of  the  Wave  Amplitude  for  Areas  Near  the  Source  and 
Near  the  Front  of  the  Wave  Train 

It  has  been  demonstrated  in  a paper  by  Whalin  (Refi  54,  60)  that  the 
asymptotic  solution  obtained  through  integration  by  the  method  of  stationary 
phase  is  valid  within  5%  provided  the  distance  r where  the  amplitude  is 
computed  is  approximately  greater  than  5 times  the  radius  of  the  initial 
surface  deformation.  If  one  desires  to  predict  the  wave  train  closer  to 
the  source  of  disturbance,  an  alternate  integration  technique  is  necessary. 

The  method  described  in  this  section  cons  ists  of  fitting  a fourth  degree 
polynomial  to  each  half  oscillation  of  the  integrand  and  to  then  perform 
the  integration  directly. 

The  solution  can  be  represented  by 

co 

r)  (r,  t)  = J f ( a)  d cr 

o 

a 3 °° 

= j f (a)  do  + J f ( a ) d a + jf(a)da 

o a 3 

where  a = a small  number  such  that  the  oscillatory  function  of  the 

integrand  may  be  approximated  by  asymptotic  expansion  near  0 = 0 

3 = a large  value  of  a such  that  the  latter  integral  is  negligible 

and  is  a zero  of  the  integrand 

ft  >t)  is  an  oscillatory  function  of  o depending  upon  the  nature 
disturbance.  Havi  ng  chosen  the  proper  limits 
a and  3,  the  first  integral  is  integrated  directly  and  the  latter 

is  negligible. 

Let  a;  be  the  ordered  zeros  of  the  function  £(»)  in  the  interval 
a < o i * 3,  i = 1,  2,  ..,  N and  0 N = 3.  then 
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f ( o)  = K {a.'  a4  + B.'  a,  + C!  a,  + D.’  a + E.’  } 
t 1 li  l Z l i J 


a £ a £ a^,  i = 0 


ai_1  < a ^ i = 2,  3,  . . . , N 


where  K is  a constant  depending  on  the  initial  conditions 

Upon  substitution  of  the  above  expression  for  f(o)  and  subsequent 
term  by  term  integration,  the  wave  amplitude  becomes 

n(r.t)  „ K +b;  <a  7 « - a7'2,  ♦ <:;<<>*«-  »5/2> 


D1  3/2  3/2. 


1/2  1/2, 


+ -y-  (a  j""  - a ' ) +E{  (a  - a"“) 

+ [t  (ai/2  - ai!f)+T  (ai7/2-  ai7-/2)  +x(ai5/2-  ai5-i2) 

i=2 

4<03/2.„3/2)+e,(0./2.01/M} 


where 


= <°i  + ai.i>Ai  +Bi 

* °i  Vl  Ai  ’ (°i  Bi  +Ci 

= °i  °i.l  Bi-<0i  + ai-l)Ci 


* °i  a,-l  Ci 


( °i  ■ °i.iy 


[|  i - 1 +| , RP=*)] 


,3  o . + a. 


A 


r 


B.  = 


5 0.  +3  o:  a +30  . 


1 ( °.  - O.  )' 

' 1 1-1' 


.(—6 


4 ) 


/ \ / + a-,-  n /3  o + 5 o . ,3a.  + a. 

. ( a.  + a . } f j -1 ill ) + ( i f 1 

^ i i-l'  \2  J \ 6 / 1 V 4 J J 


‘ ' (a  - a ,)4  5 £ V 4 J 

' i i-l' 


( o.  + 3 a . ) (3  a.  + a.  ) a.  + a. 

i i-l'  i i-l  , f i i-l > 

£ { — j — ; 


2 ( a . + 3 a . ,)(a.+  a..  .)  .3  a + a lS 

' i i-l  i i-l  , f i i-M 

+ 3 f 'v 4 J 


B,  - a . A, 

1 i 1 


C1  * a 1 B1 


D1  - a i ci 


- a i Di 


1 r<f5‘f7>  (f6-f7h 

A1  " (a  - b)  L (a  - c)  (b  - c)  J 


1 Ra  + b + c)  ,r  (b  +c  + d)  lf  . ( {] 

B1  = (TTb)  L (b  -.c)  (f6  - f7  • (a  - c)  a f7>J 


r _ (be  + b d +cd)  f _ (ac  + a d ted)  f + (*b  +&d  tb4)  f 

C1  _ (a-b)  (a-c)  f5  (a  - b)  (b  - c)  6 (a  - c)  b - c)  7 


TIME,  sec 

Fig.  2.  Wave  amplitude  as  a function  of  time  10  km  from  the  source  for  a parabolic  im- 
pulse. II  = 1.4  km,  h = =5  km,  U — 5 X 10r  dynes  sec/cm“.  Method  of  stationary  phase  (solid 
line),  integration  method  of  this  paper  (dashed  line). 


TIME,  sec 


I imc.,  act- 

Fig.  3.  Wave  amplitude  as  a function  of  time  5 km  from  the  source  for  a parabolic  im- 
pulse. R = 1.4  km,  h = 5 km,  /„  = 5 X 10’  dynes  sec/cm’.  Method  ol  stationary  phase  (solid 
line),  integration  method  of  this  paper  (dashed  line). 
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a 


--(b-a2d)  (c-a3d)-)  r(a  - d)  (c  - a3  d) 

D1  = a L b (b  - d)  + c (c  - d)  j + b L a (a  - d)  ' c (c  - d) 


r(b  - d)  (a  - d)-, 

+ c L b (b  - d)  “ a (a  - d)  _ 


<7  ^ + 3 a 

£ =£‘-f< 
5 a - d 

: o ^ +a 

f = f2  ' f4 
6 b - d 

3 a +a 
4 

l = f3  ''*4 
*7  c - d 

_ 

0 - a 

“1  * £l/£4 

0 . + 3 a 

- 4 f (—4 — ) 

a2  = f2/f4 

3 (°  l +a) 

a3  = f3^f4 

a , + a 
2 > 

a i - a 


3 c,  + a 
- 4 f (— j ) 


Oj  - a 


Examples  of  this  integration  method  are  shown  in  Figures  1 to  3 for 
the  case  of  a parabolic  impulse  and  further  details  can  be  found  in 
Whalin  (Ref.  60  ).  The  method  is  relatively  easy  to  use  and  the  accuracy 
of  the  integration  procedure  is  well  established. 


2.5  Ener 


in  Explosion  Waves 


One  of  the  important  considerations  is  the  efficiency  of  an  explosion 
as  a wave  producing  mechanism.  Many  experimental  tests  have  indicated 
that  approximately  1-3%  of  the  available  energy  appears  in  the  form  of 
surface  waves  for  HE  detonations.  The  percentage  is  considerably  more 
uncertain  in  the  case  of  nuclear  explosions.  This  section  gives  the 
expressions  for  computing  the  potential  energy  of  three  surface  deforma- 
tions assumed  for  purposes  of  predicting  the  water  waves.  A derivation 
due  to  Sakurai  is  included  for  approximately  relating  this  expression  to 
field  measurements  of  the  envelope  maximum  at  various  stations,  a dis- 
tance r from  GZ.  Further  theoretical  expressions  and  curves  are  given 
showing  the  relationship  between  the  potential  and  kinetic  energy  appearing 
in  the  wave  train  at  any  time  t after  the  detonation. 
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R 

2 r r 
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The  potential  energy  of  the  deformation  with  lip  represented  by 
r\o{ro)  = 


C nr4  4 r 2 
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A third  case  of  interest  is  an  initial  deformation  in  the  form  of  a para- 
bolic cavity  surrounded  by  a lip  which  is  an  extension  of  the  cavity.  This 
initail  deformation  is  represented  by 


% (ro> 


R 


r ^ R 
o 


r > R 
o 
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and 


,R 


E = P gr\. 


4 r 
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r dr 
o o 
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= 17  P g^0 


4R2  . 4R2  ' 

— + -r~J 


2 D2 

^ P g ’Hr,  R 


It  should  be  noted  that  the  above  surface  deformation  conserves  the  mass 
of  the  fluid  domain,  i.  e. , the  amount  of  water  above  the  undisturbed  free 
surface  is  equal  to  the  volume  of  the  water  cavity. 

The  above  deformations  are  considered  since  they  are  commonly  used 
in  performing  predictions  of  the  surface  waves  generated  by  underwater 
explosions.  Many  other  shapes  and  conditions  could  be  used;  however, 
the  present  state  of  the  art  has  indicated  that  these  initial  conditions 
yield  the  most  reliable  prediction. 

2.5.2  Kinetic  and  Potential  Energy  Distribution 

At  time  zero,  the  mathematical  models  used  for  performing  predictions 
assume  that  all  the  energy  partitioned  to  the  water  is  in  the  form  of  poten 
tial  energy.  One  can  study  the  partitioning  between  potential  and  kinetic 
energy  as  a function  of  time  after  detonation  from  the  theoretical  solution 
for  the  velocity  potential.  However,  at  large  times  after  the  explosion, 
the  total  energy  is  equally  divided  between  kinetic  and  potential  energy 

which  is  easily  demonstrated  as  follows. 
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where 


A precise  distribution  of  the  potential  and  kinetic  energy  as  a function  of 
time  can  be  evaluated  from  the  expression  for  the  velocity  potential.  The 
limiting  value  of  an  equal  partitioning  of  the  kinetic  and  potential  energy 
will  be  reached  relatively  quickly. 


2.5.3  Relationship  Between  the  Wave  Envelope  and  Initial  Conditions 

This  section  discusses  an  approximate  method  of  estimating  the  explosive 
energy  converted  to  water  waves  from  field  measurements  of  the  wave 
amplitude  as  a function  of  time  at  any  location  r from  SZ  in  water  of 
constant  depth.  The  method  presented  was  first  developed  by  Sakurai 
(Ref.  51)  and  greatly  facilitates  an  estimate  of  wave  energy. 

The  theoretical  solution  for  the  wave  amplitude  which  yields  very  good 
agreement  with  field  measurements  provided  the  cavity  parameters  qQ 
and  R are  well  chosen  is  given  by 


T)  (r,  t)  = 


l J (SiL)  / 0 ?M~~  cos  (or  - Jo  tanh  at)  (para- 


2 ' h ' v -cp'  (a) 


bolic  cavity) 


, . 4rioh  T ,/r  ii  ) /cola)  (nr  - ./rrtanh  gtl  (era- 

(r,  t)  = y J4  W J h V - Cp'  (0  ) 


ter  with  lip) 


. R)  /CTcp<oT  cos  (or  - Jo  tanh oT  (parabolic 
T)  ( r , t)  = J3  (OR)  V-cp’(a) 


depression  and  cylindrical  elevation) 
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It  is  very  easy  to  perform  measurements  of  the  first  envelope  amplitude 
from  wave  records,  which  forms  the  basis  of  this  method.  The  product 
of  the  first  envelope  amplitude  and  the  distance  r of  the  measurement 
station  is  given  by 


[parabolic  case] 


T)  r 
'r 


m m 


[lip  case] 


■H  r = n R — — (a„R)  ./°m  (parabola  and  cylinder) 

e o <7  3 v 7“7 r 

e m 

where 

a = the  value  of  a at  the  envelope  maximum 
m 

The  validity  of  this  method  is  subject  to  the  following  assumptions: 

1.  The  measurement  station  is  a large  distance  r from  the  explo- 
sion so  that  the  above  equations  are  valid. 

2.  The  value  of  o at  the  envelope  maximum  is  sufficiently  large  so 
the  asymptotic  expression  may  be  used  for  the  functions  cp(o) 
and  cp'  (a)  . 

3.  The  wave  train  is  represented  relatively  accurately  by  the  above 
equations  provided  proper  deformation  parameters  qQ  and  R 
are  chosen. 
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Parabolic : 


t / amR  \ 

rj  r = 2VT  rj  R — — R-h  ■ - 

e 'o  , a R N 

( m \ 


V3  amR 


Grater  witn  Lip: 


n*r = ^ drfr 

( N m 


Cylinder  and  Parabola:  t]  r = ^2~  nQR  (cr^R) 


Parabolic  case:  The  first  maximimum  of  rier  is  given  by  the  first 

maximum  of  J^(x)/x  , given  by 


a R 


= 2.2999 


T ° R 

a R 
(-2-) 


= 0.17996 


Therefore 


ri  r & 0.  509  rioR 


nQR  « 1.964  r\&R 


From  the  previous  section,  the  potential  energy  of  a parabolic  cavity 


is  given  by 
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Upon  substitution 


3.  86  tt  p 


2 2 

ti  r 

a v 


E«  129  n 2 r2  ft.  lbs. 
'max 


Lip  Case: 


In  this  case  the  maximum  of  J^(x)/x  is  0.0787,  resulting  in 


•n  r = 0. 771  ti  R 
'max  o 


noR  = 1.299  r!maxr 


The  total  energy  of  the  initial  deformation  is  given  by 


E « 66.  2 n 2 r2  ft*  lbs* 

'max 


Cylinder  and  Parabola  Case: 


The  maximum  of  the  wave  envelope  is  merely  at  the  maximum  of 


J3(x)  which  gives 


ti  r = 0.  614  n R 
'max  o 


%R  = 1-63Vaxr 


2 . 6 Initial  Conditions 


It  has  been  demonstrated  in  the  previous  section  that  one  has  a 
wide  variety  of  initial  conditions  in  order  to  simulate  the  wave  generation 
mechanism  for  an  explosion.  Intuitively,  for  the  case  of  an  air  explosion, 
a time  dependent  pressure,  P (rQ,  t),  would  seem  to  be  the  logical  choice 
of  an  initial  condition.  Unfortunately,  to  the  knowledge  of  the  authors,  there 
have  been  no  theoretical  developments  other  than  those  presented  in  the 
previous  section.  It  has  been  demonstrated  (communicated  through  private 
discussion  with  Drs.  Van  Dorn  and  Olsen)  that  the  integration  of  pressure 
measurements  from  field  tests  to  give  an  initial  impulse  and  the  subsequent 


use  of  an  initial  impulse  for  predicting  waves  results  in  an  erroneous  prediction 
by  more  than  an  order  of  magnitude.  The  use  of  a time  dependent  pressure 
is  the  next  logical  step  toward  a realistic  model  and  would  comprise  both 
an  interesting  and  valuable  research  task. However , tn  this  report  we  are 
interested  in  underwater  explosions,  primarily  because  the  energy  converted  to 
water  waves  is  significantly  larger  than  for  air  explosions.  It  has  been  well 
established  that  the  use  of  an  impulse  or  a surface  deformation  as  an  initial 
condition  results  in  a wave  train  which  is  qualitatively  similar  to  those 
measured  from  explosions.  After  considerable  trial  and  error,  it  was 
further  demonstrated  that  a parabolic  deformation  or  impulse  of  the 


i 


ii  , uJ 


resulted  in  the  best  comparison  with  experiment  provided  the  two  parameters 
T\  » R or  ^0»  R were  properly  selected.  The  only  objection  to  this  form 
of  initial  condition  is  the  fact  that  an  initial  deformation  consisting  of  a 
net  addition  or  removal  of  fluid  from  the  domain  results  in  the  prediction  of 
a discontinuity  in  the  water  surface  at  the  front  of  the  wave  train.  There- 
fore any  prediction  of  the  leading  wave  will  be  in  error,  and  under  certain 
conditions,  the  leading  wave  is  of  considerable  interest  due  to  its  long 
period  and  large  shoaling  coefficient. 

Whalin  (Ref.  57  ) developed  a model  consisting  of  a water  crater 
with  a lip  where  the  volume  of  fluid  in  the  domain  is  conserved,  thus 
eliminating  the  discontinuity  at  the  front  of  the  wave  train.  This  deformation 
is  given  by 

r £ ./T  R 
r > J~2>  R 


1 'o  (ro}  fro  L 0 „4 


4 r 


+ 1 


3 R 3 R 


and  is  illustrated  in  the  following  figure. 


7 7 7 7 7 7 7 7 7 7 7 7 7 *7  7 


Figure  4 Initial  Deformation  with  a Lip 
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Recent  preliminary  investigations  indicate  that  perhaps  the  com- 
bination of  a parabolic  deformation  and  cylindrical  elevation  gives  a better 
prediction  of  the  entire  wave  train  than  either  of  the  previous  two  initial 
conditions.  This  deform  ion  is  given  by 


Figure  5 Cylindrical  Elevation  plus  Parabolic  Depression 

Therefore,  through  a considerable  amount  of  investigation,  the  initial 
conditions  which  are  commonly  used  for  predicting  water  waves  generated 
by  underwater  explosions  consist  of  one  of  the  stationary  surface  deforma- 
tions above. 

Several  investigations  need  to  be  conducted  to  take  lull  advantage 
of  the  presently  available  theoretical  developments.  These  include: 
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1.  The  accuracy  of  the  above  initial  conditions  for  predicting  wave 
amplitudes  near  the  front  and  end  of  the  wave  envelope.  It  is 
well  established  that  periods  and  maximum  amplitudes  can  be 
predicted  accurately. 

2.  The  applicability  of  time  dependent  deformations  or  an  initial 
velocity  field  for  developing  a better  model  of  the  initial  con- 
ditions. 

3.  The  accuracy  of  predictions  close  to  the  explosion,  yet  away 
from  the  base  surge  area. 

All  of  the  above  investigations  require  a close  comparison  between 
theory  and  field  tests.  Even  more  important,  field  tests  are  a necessity 
in  evaluating  the  initial  conditions  r]  and  R as  discussed  in  the  fol- 
lowing section. 


The  previous  theoretical  developments  are  all  based  on  the  linear 
theory  of  surface  waves  in  water  of  constant  depth.  Obviously,  the  hydro- 
dynamic  phenomena  in  the  immediate  vicinity  of  an  explosion  constitute  a 
highly  nonlinear  problem.  Initially,  it  seems  absurd  to  expect  a linear 
theory  assuming  a mere  deformation  of  the  water  surface  to  yield  reliable 
predictions  of  the  water  waves  produced.  However,  through  comparisons  of 
the  theoretical  results  with  field  tests,  it  was  well  established  that  the 
linear  theory  produced  qualitatively  satisfactory  results.  The  only  problem 
remaining  was  to  relate  the  parameters  r]  and  R of  the  initial  deformation 
to  charge  weight,  W,  submergence  depth,  z,  and  water  depth,  h.  In  the 
case  of  explosions  in  deep  water,  defined  as  water  sufficiently  deep  such 
that  the  bottom  has  a negligible  effect  on  the  generation  of  the  surface  waves, 
the  only  parameters  of  consequence  are  W and  z. 

The  present  state  of  the  art,  which  seems  to  be  quite  satisfactory 
for  logistic  purposes,  could  not  have  been  achieved  without  the  performance 
of  this  particular  experimental  test  series.  The  range  of  charge  weightsand 
submergence  depths  used  in  this  test  series  was  sufficient  to  result  in  a 
reliable  empirical  evaluation  of  the  parameters  'H0(z»  W)  and  R(z,  W). 

This  initial  evaluation  of  these  parameters  was  performed  by  Whalin  (Ref.- 57  ) 
without  all  the  data  from  this  test  series.  These  parameters  are  reevaluated 
(Ref.  31  ) in  Section  4 of  this  report  with  the  inclusion  of  all  data  from  the 
program.  The  method  of  evaluation  is  presented  in  the  following  paragraphs. 

A dimensional  analysis  of  the  problem  was  attempted, 
however,  this  procedure  failed  due  to  the  lack  of  any  consistent  plot  of  the 
data  in  terms  of  dimensionless  parameters.  This  was  probably  due  to 


many  factors:  the  fact  that  gravity  could  not  be  scaled,  a change  in  efficiency 

of  detonation  with  increasing  charge  weight,  the  speed  of  detonation  could  not 
be  scaled,  etc. 


Approaching  the  problem  from  an  empirical  point  of  view,  the 
following  steps  were  undertaken: 

1.  A confirmation  of  the  1/r  amplitude  decay  of  the  envelope  with 
distance.  This  has  been  verified  (Refs.  57,  50,  51  ). 

2.  An  empirical  formula  for  the  period,  wave  number,  and  wave 
length  associated  with  the  envelope  maximum  (all  are  related  to 
each  other)  as  a function  of  z and  W . 

3.  An  empirical  relationship  for  the  product  of  the  envelope 
amplitude  and  distance  from  SZ  (Tlmax  r which  was  shown  constant 
by  1.  above)  as  a function  of  z and  W for  any  one  detonation.  . 

The  relationships  above  can  be  represented  by 


k = f,  (z,  W) 
max  1 


p r = f?  (z,  W) 

'max  2 

order  to  evaluate  n and  R from  these  relationships,  one  must  refer  to  the 


asymptotic  solution  for  the  case  of  an  initial  deformation 


n (r,  t)  - 14^  cos  r - J 0 tanh  0 *> 


where  the  envelope  is  represented  by 


"e  (r,  t) 


Va) 


fa  cp  (o  F 
- cp ' ( a) 
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Note  that  nQ  (a)  contains  both  parameters  T]o  and  R and 


(J) 


-g—  (a  R)  parabolic  deformation 


T]o  (a)  = — 2~  (a/3  oR)  deformation  with  a lip 


% (0)  = 


— — (a  R)  cylinder  and  parabola' 


We  now  proceed  with  the  evaluation  of  R.  Noting  that  g = kh  or 

a = k h = hf.  (z,  W),  the  value  of  a at  the  envelope  maximum, 
max  max  1 

d(r|  r) 

0 , is  the  solution  of  — -j— = 0 . Therefore,  R is  the  solution  of 


r d(r!er)  i 

L d a ja  = a = hf.  (z,  W) 

u u max  1 


0,  shown  below  for  the  specific 


cases  aDove. 


d / 1 T , / OQ  (g)  v 1 

do“  1 2 J2  (aR)  V -cp'(a)  / J a = hf^z,  W) 


= 0 parabolic  case 


}"  d / 1 

L d a C 2 

a 


h }]0  = h£1(z.W)  = li‘> 


L®-  {’TJ3toR)  J^&)  i J a = hfj  U. 


gffl(g)  11  _ =0  parabola  and  cylinder 


The  solution  for  n is  now  a simple  matter,  upon  substitution 

' n 


. . Jk 


parabolic  case 


T| 


O 


f2(z'W)  T a2 
2 L J2  (a  R) 


/-cp'(p)  1 

V acp(a)  J a = hfL  (z,  W) 


% 


f2  (z.  W) 
4 


r SL HEM  1 

L J4cyi"aR)  ' acp(a)  Joshf1(*iW) 


lip  case 


O I -CD1  (c)  1 

no=f2(z,w)  LRT3(oR)  Ja  = hf1(z.W)  cylinder  and  parabola 

The  functions  and  are  evaluated  in  Section  3 of  this  report 
along  with  a table  of  n and  R as  a function  of  z and  W . 

One  must  realize  at  this  point  that  the  values  of  r]  and  R obtained 
in  this  manner  will  not  necessarily  have  any  resemblance  to  the  actual  water 
crater  formed  by  the  explosion.  A theoretical  solution  has  been  developed 
which  adequately  describes  the  qualitative  features  of  the  wave  train.  The 
experimental  data  are  then  used  to  evaluate  the  initial  conditions  of  this 
model  as  a function  of  z and  W so  that  the  predictions  will  result  in 
empirically  derived  values  of  the  wave  number  and  amplitude  at  the  en- 
velope maximum. 

It  is  readily  apparent  that  reliable  predictions  could  not  be  performed 
without  both  the  theoretical  analysis  of  problem  and  the  performance  of 
well  planned  and  executed  field  tests. 
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